Valence-bond crystal in a {111} slice of the pyrochlore antiferromagnet 
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We investigate theoretically the ordering efltect of quantum spin fluctuations in a Heisenberg anti- 
ferromagnet on a two-dimensional network of corner sharing tetrahedra. This network is obtained as 
a {111} slice of the highly frustrated pyrochlore lattice, from which it inherits the equivalence of all 
three pairs of opposite bonds of each tetrahedron. The lowest-order (in 1/S) quantum corrections 
partially lift the huge degeneracy of the classical ground state and select an ensemble of states with 
long-range valence-bond order. 



Magnets with strong geometrical frustration have be- 
come a focus of experimental studies and inspired a grow- 
ing number of theoretical inquiries into their unusual 
properties^ For classical spins, strong frustration pre- 
cludes simple Neel order and creates a vast ground state 
degeneracy. This degeneracy, in turn, renders the mag- 
net susceptible to nominally small perturbations (dipolar 
interactions, anisotropics, quantum effects etc.). 

An antiferromagnet on the pyrochlore lattice — a net- 
work of corner-sharing tetrahedra — is perhaps the ulti- 
mate example of strong frustration. Monte Carlo simu- 
lations with classical Heisenberg spins show that thermal 
fluctuations alone are ineffective in restoring magnetic 
order: the magnet remains in a paramagnetic but cor- 
related state at a temperature as low as 10^'' JS*^. The 
effect of quantum fluctuations on the pyrochlore antifer- 
romagnet has not yet been determined. 

Previous studies have focused on the case of maximal 
quantum fluctuations, S = 1/2. Early on, a ground 
state in the form of a valence-bond crystal (VBC) was 
proposed,^ which preserves the spin-rotational 0(3) sym- 
metry but breaks a number of discrete lattice symmetries, 
including the equivalence of the three pairs of opposite 
bonds of a tetrahedron. The relevant local order param- 
eter is constructed from the expectation values (S^ • Sj). 
This result was obtained, in the absence (to this day) 
of a demonstrably reliable analytical treatment, via an 
uncontrolled approximation; several methods have since 
yielded similar results.^' 

An alternative route of studying quantum effects be- 
gins in the classical limit 5 — > oo. At large S, quantum 
fluctuations are weak and can be evaluated in the frame- 
work of a perturbation theory, wherein physical quanti- 
ties are expanded in powers of 1 /S. The large-S" approach 
has been previously applied to a variety of strongly frus- 
trated magnetsAiiSiS 

In this paper, we report an application of the large-S" 
method to a {111} slice of the pyrochlore lattice shown 
in Fig. n^a). The resulting lattice can be described as 
a kagome plane flanked by two triangular lattices, so 
that every triangle is promoted to a tetrahedron. (Note 
that this is different from the pyrochlore slab geometry 
of SCGOjlS where a triangular lattice is sandwiched be- 



tween two kagome planes.) 

From the symmetry point of view, our {111} slice has 
the following features important in connection with the 
pyrochlore in 3D. All its units are equivalent tetrahe- 
dra and it thus permits coUinear classical ground states 
(unlike SCGO). Further, unlike the checkerboard version 
of the 2D pyrochlore,^ it retains equivalence between all 
three pairs of opposing bonds of each tetrahedron, a sym- 
metry spontaneously broken in the 5* = 1/2 pyrochlore 
VBC mentioned above 4 

For this lattice, we have found and characterized the 
subset of classical ground states favored by quantum fluc- 
tuations. The problem of minimizing the zero-point en- 
ergy is equivalent to a discrete gauge-like theory, which 
by its very nature has many degenerate ground states. 
A natural flrst guess is that the ground states are those 
with uniform flux; this ensemble leads to a semiclassical 
valence bond liquid. However, our numerical minimisa- 
tions reveal vacua with a pattern of gauge fluxes violat- 
ing translational symmetry; this is reflected in a spatial 
modulation of nearest-neighbor spin correlations (Si - Si). 
Along with the absence of Neel order, this implies that 
the spin system is a valence-bond crystal. 

We study the Heisenberg spin Hamiltonian 

H = S,-S, = (J/2) J2 |La|'+const = 0{S^), (1) 

(ij) " 

where J > and the sum is taken over all bonds {ij) 
shown in Fig.^a). In the alternative representation, the 
Greek index a enumerates tetrahedra; Lq is the sum of 
the four spins residing on tetrahedron a. The energy 
is minimized when Lq, = for every tetrahedron a. Clas- 
sically, this condition does not determine a unique ground 
state. (This is the hallmark of frustration.) Rather, it de- 
fines an enormous manifold of classical vacua. By explicit 
construction, or Maxwellian counting)^ the dimension of 
the manifold can be seen to equal twice the number of 
tetrahedra. Thermal order by disorder is absent. 

Quantum fluctuations produce virtual excitations out 
of the ground state, thereby lowering the vacuum energy. 
For large S, the quantum correction can be expanded in 
powers of 1/S' by expressing the spin operators in terms 
of Holstein-Primakoff bosons. (The virtual excitations 
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FIG. 1: (a) A generic coUinear ground state of Eq. Q on 
the pyrochlore slice. Open (filled) dots correspond to Si — 
+1 (—1). (b) Kagome sites live on links of a honeycomb 
lattice. The Z2 flux is —I through every hexagon, (c) and (d) 
Resulting dimer coverings of the two triangular sublattices. 
The dimers cross frustrated bonds {si = sj). 



are magnon pairs produced by the operators feffoj in the 
bosonized Hamiltonian.) The lowest-order correction is 
the zero-point energy of harmonic spin waves, 



'2NJS + M^a\/'2 = 0{S), 



(2) 



where {ua} is the set of spin- wave frequencies in a given 
classical ground state and N is the number of tetrahe- 
dra. The first quantum correction |(2Jl depends implicitly 
(through the magnon spectrum) on the spin configura- 
tion and thus can be viewed as an effective spin potential 
breaking the degeneracy of the ground-state manifold. 

Typically, quantum fluctuations select coUinear ground 
stateSfii provided that such vacua exist at the classical 
level. Although we are unaware of any general proof of 
this proposition, a compelling argument in its favor can 
be stated nonetheless. Given that the reduction of the 
vacuum energy is due to generation of virtual magnon 
pairs and that magnons are transverse excitations, the 
zero-point energy should be lower when more adjacent 
spins share both of their transverse directions, i.e. when 
spins are coUinear. fn certain cases this heuristic argu- 
ment can be made preciseiSiiS We assume that it applies 
in this case as well and minimize the zero-point energy in 
the set of coUinear classical vacua. As a quick check, we 
have computed zero-point energies of a few simple classi- 
cal vacua with the smallest magnetic unit cell, and found 
them to be -1.2126JS', -1.1612JS', and -1.1436JS' for 
the coUinear, "perpendicular," and "anticoUinear" states, 
respectively (notation of Ref . ITll) . 

We parametrize coUinear states by means of Ising vari- 
ables Si = ±1 such that Si/ S = SiU, where n is a unit 



vector. The constraint of zero total spin on every tetrahe- 
dron carries over to the Ising variables: so + si +S2 + S3 = 
0. We can eliminate the variable sq on the out-of-plane 
site at the expense of introducing another constraint: 



si + S2 + S3 = ±1 on every triangle. 



(3) 



It follows from Eq. @ that coUinear classical vacua of 
the pyrochlore slice can be mapped onto the ground 
states of a kagome Ising antifcrromagnct. Redundancy 
of the out-of-plane Ising variables allows us to speak in- 
terchangeably of triangles and tetrahedra when we deal 
with ground-state configurations. 

The zero-point energy |(2J) is an implicit function of 
the Ising variables {si}. As outlined in Refs. |0 and 
1^ the spin-wave spectrum can be obtained by lineariz- 
ing the classical equations of motion for the spins, writ- 
ten in terms of the L variables: KLa — X]/3 ^ 
Here Sap denotes the spin shared by tetrahedra a and 
p. (All of these spins reside in the kagome plane.) 
For a coUinear spin state polarized along n = (0,0,1), 
transverse fluctuations of the total spin on tetrahedron 
a can be parametrized by a complex number <Ta as 
Jja — (Rsctq, ImCTc, 0). The magnon frequencies are then 
determined by solving the eigenvalue problem 



JS* V Sq/3 cr/3 = huaa- 



(4) 



The zero-point energy can now be minimized numer- 
ically, e.g. by using simulated annealing. However, it 
turns out that the states minimizing Eq. |2Jl retain a 
large residual degeneracy. Hence, finding one such state 
does not yet solve the problem. We must first understand 
the nature of the residual degeneracy. This is best done 
by discussing a hidden symmetry responsible for it. 

A gauge-like symmetry. Observe that the triangles 
{a} occupy the sites of a honeycomb lattice dual to the 
kagome [Fig.^b)]. The in-plane sites {af3} live on links 
of the dual lattice. It was noted by Henley^^ that Eq. Q 
reveals a special gauge-like symmetry of the zero-point 
energy |(2Jl. Namely, a transformation 



Sap 



= AqSq^Ao"^ (5) 



leaves Eq. Q invariant. (Conservation of Ising spin 
length \sai3\ requires that A^ = ±1.) Therefore, two 
spin configurations {saf}} and {s^^} related by the gauge 
symmetry (jS)) have identical magnon spectra and, conse- 
quently, equal zero-point energies. We therefore expect a 
large residual degeneracy of the ground state even after 
the inclusion of the first quantum correction (0). 

A word of caution. Transformation jSJl is not exactly 
a gauge symmetry, even though it strongly resembles 
one. The transformed state must respect the ground- 
state constraint Q. The importance of this restriction 
is evident when one ponders the meaning of the trans- 
formation: if Aq = —1, we Uip all Ising spins on trian- 
gle a. Doing so may violate the ground-state constraint 
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(O on adjacent triangles, where only one spin is flipped. 
This reduces the number of possible gauge transforma- 
tions but — crucially — does not eliminate them entirely. 

The uncovering of a hidden symmetry greatly simpli- 
fies our task: as soon as we find a single ground state 
minimizing the zero-point energy (j^J, we can generate 
all of the others by applying the gauge transformations 
But the gauge- like symmetry is useful in more ways 
than one. It also puts restrictions on the way in which 
the zero-point energy depends on the Ising spin variables. 

From the viewpoint of the Z2 symmetry (UJ , the Ising 
spins Sq/3 are gauge variables. In a true gauge theory, the 
energy can only depend on their gauge-invariant com- 
binations, in this instance the Zi fluxes through closed 
loops of various lengths I on the hexagonal lattice, 
<\>it = Saia2 SQ2a3 ■ • ■ Sqjqi • The flux through any such 
loop can be factored into a product of fluxes piercing 
elementary plaquettes of the dual lattice, the hexagons 
7 = Fg [Fig. ^b)]. Therefore, the zero-point energy will 
be a function of Z2 fluxes: 

i?i (6) 

While ours is not precisely a gauge theory, Eq. © is 
the right guess: the zero-point energy is entirely deter- 
mined by the values of the Z2 fluxes. To prove this, we 
demonstrate that any two coUinear ground states with 
the same set of flux variables have equal zero-point 
energies (0. To accomplish that, we show that the two 
states are related by a gauge- like transformation (jSJ. 

Proof. Let the Ising link variables be {s^^} and {s^^} 
in the two states. Deflne new link variables as their ratio, 
tai3 = s'^/s/SajS = ±1- The variables tajs are pure gauge: 
their product along any closed loop equals There- 
fore we can unambiguously define a new variable on 
sites of the dual lattice in such a way that Aq, = tapAp 
everywhere. Hence s'^^ = AqSo/jA^^. By Eq. the 
zero-point energies of the two states are equal. Q.E.D. 

The precise form of the flux potential is not fixed by 
Eq. © and can be rather complicated. One might hope 
that a series expansion could successfully mimic it: 

7 7,7' 

(7) 

If the linear term dominates, the ground state has uni- 
form fiux (pry = — sign(ai) on all elementary plaquettes. 
On the checkerboard, cj)^ — +1 in the ground statei^ 

Extending the same simplistic approach to the py- 
rochlore slice we can take a guess that the ground states 
will either have no flux (0^ = +1) or have the maximal 
flux TT (0^ = —1). While neither of these is the correct 
answer, they are interesting in their own right. 

The TT-flux states. In these ground states, the prod- 
uct of six Ising spins around any hexagon (of the original 
kagome lattice) equals —1 [Fig. QJa)]. The tt flux on the 
hexagon is related to the number of frustrated bonds. 
(A bond is frustrated if it connects Ising spins Si = Sj.) 



More precisely, each hexagon has or 2 frustrated bonds 
belonging to the triangles pointing down [Fig.Q^b)]. This 
fact allows us to map any 7r-flux state onto a dimer cov- 
ering. If a hexagon has 2 such frustrated bonds, draw 
a dimer connecting the centers of the two triangles con- 
tributing the frustrated bonds [Fig. QIc)] . Because every 
triangle has exactly one frustrated bond, it shares a dimer 
with one other triangle. The same construction applies 
to the sublattice of triangles pointing up [Fig.^d)]. We 
note that gauge transformations on the down triangles 
are sufficient to take any of the dimer states on the up 
triangles into any other without changing the dimer state 
on the down triangles - and vice versa. Thus every 7r-fiux 
ground state maps onto a pair of completely independent 
dimer coverings of the triangular lattices. 

Now one can use the well-known properties of the 
dimer model on the triangular latticei^ to characterize 
the ensemble of the 7r-flux states. By construction, the 
number of such states equals the number of pairs of dimer 
coverings. The latter grows exponentially with the num- 
ber of sites, so that there is a flnite entropy of 0.4286 per 
tetrahedron. A disordered state of classical dimers on a 
triangular lattice translates into the absence of valence- 
bond order in our spin model. Spin correlations are also 
short-ranged. On a lattice with periodic boundary con- 
ditions (on a torus), these ground states fall into four 
disconnected topological classes distinguished by the val- 
ues of the global Z2 fluxes piercing the torus. 

In short, the ensemble of 7r-flux states describes a clas- 
sical valence-bond liquid. 

The actual ground states. We have performed a nu- 
merical minimization of the zero-point energy ^ in the 
space of collinear vacua Q using the method of simulated 
annealing. A Monte Carlo step, conserving the classical 
energy , consisted of flipping a randomly chosen chain 
of alternating spins. The zero-point energy (01 was then 
determined numerically. 

Initial runs were done on systems containing between 
100 and 200 tetrahedra. Starting from a periodic state 
we randomized the system by flipping randomly chosen 
antiferromagnetic chains until spin correlations with the 
original state were lost. The lattice was then annealed to 
achieve a state with the lowest energy. Throughout a sim- 
ulation, we monitored Ising spins Si, valence-bond vari- 
ables SiSj, and Z2 fluxes on hexagons. The low-energy 
states obtained by the annealing showed no apparent spin 
or valence-bond order [Fig.|21[a)]. However, the Z2 fluxes 
showed a clear tendency to order: 1/4 quarter of the 
hexagons had the trivial Z2 flux -1-1 and formed a regu- 
lar pattern [Fig. |2Ia)] . We denote such states as 2 x 2 
after the size of the unit cell. 

Finite-size scaling studies (with up to 10^ tetra- 
hedra) yielded the zero-point energies of — 1.2127JS', 
-1.2462JS', and -1.2497 JS' per tetrahedron for the no- 
flux, TT-flux, and 2x2 states, respectively, with the 2x2 
states the lowest. On the basis of the energy differences, 
we estimate that quantum selection of ground states 
will be effective below the temperature scale of roughly 



4 



»:---::y-----:i»-;---::y------bv---::y----"a:---^ 




FIG. 2: (a) A generic 2x2 coUinear ground state and the cor- 
responding dimer covering (only one dual sublattice is shown). 
SoUd lines depict the decorated triangular lattice. The num- 
bers ±1 indicate the Z2 flux, (b) Expectation values of the 
valence-bond variables (siSj) are encoded in line thickness. 

10^^ JS/kB- Note that the naive energy scale JS comes 
with a numerically small prefactor. 

To obtain a dimer model appropriate for the 2x2 flux 
pattern, we start with the dimers on a triangular lattice, 
which yield the 7r-flux states [Fig. ^c)] . To change the 



flux from —1 to +1 on one quarter of the hexagons, we 
decorate the triangular lattice by modifying one quarter 
of triangles [Fig. |2Ia)]. Again, every frustrated bond is 
intersected by a dimer. 

The properties of classical dimers on the decorated 
lattice can be readily computed by using Grassmann 
variables. The number of 2 x 2 ground states is ex- 
ponentially large, with an entropy per tetrahedron of 
0.4248. The probability to find a dimer on a short link 
is p = 0.4004 > 1/3. In spin language this means that 
(siSj) =p(-hl) + (l-p)(-l) = -0.1992 > -1/3 on links 
surrounding the hexagons marked +1. These bonds are 
therefore more frustrated than the lattice average. On 
the other hand, a tetrahedron placed symmetrically with 
respect to the three closest fluxes -1-1 has {siSj) = —1/3 
for all its bonds. The valence-bond averages (siSj) are 
shown schematically in Fig. |2tb) . They are nonuniform, 
therefore the ensemble of 2 x 2 states has valence-bond or- 
der that enlarges the lattice unit cell. (We note in passing 
a similarity with the valence-bond order proposed earlier 
for the three-dimensional S" = 1/2 pyrochlorej^ 3/4 of 
tetrahedra have unequal valence-bond correlations; the 
other 1/4 remain symmetric.) 

Spin correlations (siSj) between more distant neigh- 
bors decay exponentially with distance. The correlation 
length is 1.15 (the unit of length is the side of a tetrahe- 
dron). This points to the absence of Neel order. 

To summarize, we find that, to leading order in l/S, 
quantum fluctuations select a coUinear ensemble of de- 
generate ground states with additional bond order and no 
spin order: we have a valence-bond crystal. This degener- 
acy may be lifted in higher orders in l/S, as found previ- 
ously for the Heisenberg antiferromagnet on kagomeii^ii^ 
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